Staggering behavior of O"*" state energies in the Sp{4) pairing model 
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We explore, within the framework of an algebraic sp(4) shell model, discrete approximations to 
various derivatives of the energies of the lowest isovector-paired 0^ states of atomic nuclei in the 
40 < j4 < 100 mass range. The results show that the symplectic model can be used to successfully 
interpret fine structure effects driven by the proton-neutron (pn) and like-particle isovector pairing 
interactions as well as interactions with higher J multipolarity. A finite energy difference technique 
is used to investigate two-proton and two-neutron separation energies, observed irregularities found 
around the N — Z region, and the like-particle and pn isovector pairing gaps. A prominent stag- 
gering behavior is observed between groups of even-even and odd-odd nuclides. An oscillation, in 
addition to that associated with changes in isospin values, that tracks with alternating seniority 
quantum numbers related to the isovector pairing interaction is also found. 



I. INTRODUCTION 



The observed staggering of energy levels in atomic nu- 
clei requires a theory that goes beyond mean-field con- 
siderations 0. Staggering data contain detailed infor- 
mation about the properties of the nucleonic interaction 
and suggest the existence of high-order correlations in 
the collective dynamics. Most studies of staggering fo- 
cus on two aspects of the phenomena. There are discrete 
angular momentum dependent oscillations of physical ob- 
servables; namely, of Ml transitions in nucleil3|or of the 
energy levels themselves (e.g., in octupole JJlj, , su- 
perdeformed [1, 0, |1| , ground and 7 Q, 0, [3 bands in 
atomic nuclei, as well as in molecular rotational bands 
[Tl|). And then there are sawtooth patterns of differ- 
ent physical quantities (most commonly binding ener- 
gies) that track with changes in the number of particles 
in a system (both in nuclei |l2| and in metallic clusters 

El El). 

In nuclear structure physics, staggering behavior of 
the second type is observed when one changes in a sys- 
tematic way the usual nuclear characteristics such as 
proton (Z), neutron (N), mass (A) or isospin projec- 
tion {\Z — N\/2) numbers. Examples of these nuclear 
phenomena include odd-even mass staggering (OEMS) 
[llTlll^lTllTlll9y^j2ll|. odd-even staggering in iso- 
tope/isotone shifts [22, 123, and zigzag patterns of the 
first excited 2f state energies in even-even nuclei p^ . 
The staggering behavior of a nuclear observable is most 
easily seen when discrete derivatives of second or higher 
order in its variable (s) are considered. The aim of this 
approach is to filter out the strong mean-field (global) 
effects and in so doing reveal weaker specific features. 
In this way, for example, the OEMS, which is usually 
attributed to the nuclear pairing correlations, manifests 
itself in certain finite differences of the binding energies 
that can provide for a measure of the empirical pair- 
ing gap Likewise, various discrete approximation of 
derivatives (filters) of the binding energies can be con- 
sidered to investigate detailed properties of the nuclear 



structure [23, |26l |27|, |2a 123 130|, |31|, 133 ■ 

In this paper, we consider the binding energies of the 
0'*' ground states of eyen-A nuclei in the mass range 
40 < A < 100, and for the rest odd-odd nuclei with a 
{J^ ^ 0+) ground state, the energy of the lowest isobaric 
analog 0"*" excited state (which corresponds to the ground 
state of the even-even neighbor). We refer to these states 
as lowest isovector-paired 0+ states [s^- Our aim is to 
investigate how various, comparatively small but not in- 
significant, parts of the interaction between nucleons in- 
fluence these states when we consider higher-order dis- 
crete derivatives of their energies within the framework 
of a convenient systematics. 

The algebraic pairing model [Ul^ we exploit is based 
on a fermion realization of the syrnplectic sp(4) algebra 
which is isomorphic to so(5) [sj, yl [s^- It includes 
an isovector (isospin T = 1) pairing interaction as well 
as a diagonal (in an isospin basis) pn isoscalar (T = 0) 
part. The latter is prop ortional to a so-called T(T + 1) 
symmetry termE3 |33j. The operators of the reduction 
sp{A) D u{2) D u{l) © su{2) provide for a convenient 
and useful classification of nuclei and their corresponding 
ground and excited states. The systematics is in terms 
of the eigenvalues of these operators, namely, the total 
number n of valence nucleons and the third projection i 
of the isospin and their linear combinations. 

We have already shown in [s^ that the 5*^(4) (49| model 
leads to a good reproduction of the experimental energies 
|38|| of the lowest isovector-paired state for even-yl 
nuclei, 32 < A < 100. As pointed out [s^, although 
the T = 1 like-particle pairing energy and the T = 1 pn 
pairing energy yield An = 2 staggering patterns that are 
of opposite phases, the total isovector pairing energy has 
a smooth behavior. It is the symmetry term that makes 
an accurate theoretical prediction of the regular zigzag 
pattern of the experimental energies in isobaric sequences 
possible. As a further and more detailed investigation, we 
now consider different types of discrete derivatives of the 
Coulomb corrected |39j | energy function according to the 
Sp{A) classification and with no adjustable parameters. 

The symplectic 5^(4) scheme not only allows for a sys- 
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tematic investigation of staggering patterns in the ex- 
perimental energies of the even-A nuclei, it also offers a 
simple algebraic model for interpreting the results. More- 
over, this detailed investigation serves as a test for the 
validity and reliability of the 5^(4) model and the inter- 
actions it includes. 



II. Sp{i) CLASSIFICATION SCHEME 

We start with a brief outline of the algebraic approach 
used to interpret phenomena that have been ob- 
served experimentally and are related to the isovector 
{T — 1 pairing correlations) and isoscalar interactions in 
nuclei. The sp(4) algebra is realized in terms of creation 
'^'jma annihilation Cjma fermion operators with the 
standard anticommutation relations {cjma,Cj,j^,^,} ~ 
Sjj'Sm^m'Sa.ij' , whcrc thcsc operators create (annihilate) 
a particle of type a = ±1/2 (proton/neutron) in a state 
of total angular momentum j (half integer) with projec- 
tion m in a finite space 2Q ~ Sj(2j + 1). In addition 
to the number operator N = 7V+i + N^i and the third 
isospin projection Tq = (-/V+i — N-i)/2, the generators 
of 5*^(4) are 



T± = 



A,. 



(1) 



1 



C ■ C ■ 



(2) 



whcrc iV±i arc the valence proton (neutron) number op- 
erators, To and T± are the valence isospin operators, and 
the generators, Ag_|_]^_j^, create a proton-neutron (jyn) 
pair, a proton-proton (pp) pair or a neutron-neutron (7171) 
pair of total angular momentum J'" = 0'^ and isospin 
T ~ 1. A totally symmetric finite space is spanned by 
the basis vectors constructed as (T = l)-paired fermions 



|7^+l,r^o,n_l) = ( 



A 



Al,) " |0) , (3) 



where n+i^o,-! £^re the numbers of pairs of each kind, pp, 
pn, nn, respectively, and |0) denotes the vacuum state. 
In the like-particle pairing limit, no gives the number of 
protons (neutrons) not coupled to J ^ pp (nn) pairs 
and hence defines the usual seniority quantum number 
|40l l4l|. Vi = no- On the other hand, in the pn pairing 
limit another seniority number, 2i^o = 2n-|_i+2n_i, is rec- 
ognized that counts the particles not coupled in J = pn 
pairs. The dependence of i/q on i/i within a given nucleus 
allows one to consider only Vi in the analysis; specifically, 
for a system of n valence particles with isospin projec- 
tion i = (Z — N)/2, the fully paired states Q differ in 
their coupling mode as the seniority quantum number i^i 
(j/Q = n/2 — j/i) changes by ±2 (=f2). 



TABLE L Realizations of the w''(2) = i(''(l) su''(2) subal- 
gebras of sp(4). 



Symmetry (/i) 


(cv 


Eigenvalues 
of 




Isospin (T) 


N 


n 


T+,To,T- 


pn pairs (0) 


To 


i 


Al^N-n,Ao 


pp pairs (-I-) 


N-i 




Al„^{N+,~n),A+^ 


nn pairs ( — ) 


N+i 


N+i 


ALi,i(iV_i-n),A-i 



As a dynamical symmetry, the S'p(4) symplectic 
group describes isovector pairing correlations and isospin 
symmetry through the four different reduction chains 
5p(4) D J7^(2) D [/^(l) ® 3^(2) with n = T,0,± (Ta- 
ble |T|). The first order invariant of u''{2), cp{^'0'±> = 
{N,To,N^i}, realizes the u''(2) D sm'"(2) reduction and 
reduces the finite action space into a direct sum of unitary 
irreducible representations (irreps) of U^(2) {{n,i,N^i} 
multiplets). Within a multiplet the third projection gen- 
erator of SU^{2) (middle operator in the fourth column 
in Tabic nj further reduces the U^^{2) representation to 
a vector with fixed quantum numbers {n,i)^ or alterna- 
tively (A^+i,7V_i), to which corresponds a given nucleus 
(a cell in Table lll|) . In this way the dynamical Sp{A) sym- 
metry provides for a natural classification scheme of nu- 
clei as belonging to a single-j level or a major shell (multi 
j), which arc mapped to the algebraic multiplets. This 
classification also extends to the corresponding ground 
and excited states of the nuclei including their isovector- 
paired 0+ states. 



TABLE IL Classification scheme of even-A nuclei in the I/7/2 
shell. The shape of the table is symmetric with respect to the 
sign of i and n — 2Q. An = 2 in each i multiplet (columns), 
Ai = 1 in each n multiplet (rows), AN±i = 2 in each N^n 
multiplet (diagonals). The subsequent action of the SU^{2) 
generators (shown in brackets) constructs the constituents in 
a given SU''(2) multiplet {/i = T, 0, ±). 
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22Ti20 


210C21 


2oCa22 
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44, r 
23 V21 


22Ti22 


44c„ 
21'3C23 


2oCa24 






6 




24Cr22 


23V23 


22Ti24 


21'5C25 


2oCa26 




8 




iMn23 


48 


48 -ir 
23 V25 


22Ti26 


21'5C27 


2oCa28 


10 




igFe24 


iMn25 


24Cr26 
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23 V27 


22Ti28 




12 
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The general model Hamiltonian with Sp{A) dynamical 
symmetry, which consists of one- and two-body terms, 
can be expressed through the Sp{4) group generators. 



-GAIAo-F{AI^A+i+AI^A^i)-E'- 2^ 
^)(T2-|)-,7V, 



-C 



JV(jV-l) 



(4) 



where G, F, E, C and D arc strength parameters and 



3 



e > is the Fermi level energy. This Hamiltonian con- 
serves the number of particles and the third isospin pro- 
jection and changes the seniority quantum number vi 
by zero or ±2, the latter implies scattering of a pp pair 
and a nn pair into two pn pairs and vice versa. The 
isospin breaking Hamiltonian includes an isovector 
(T = 1) pairing interaction (G > 0, F > for at- 
traction) and a diagonal (in an isospin basis) isoscalar 
(T = 0) force, which is related to a symmetry term (E). 
Within a shell (a single-j level, I/7/2, or a major shell, 
l/(|)2p(-| i)l(?(|)), a reasonable estimate for the param- 
eters in the Hamiltonian Q is obtained in a fitting pro- 
cedure of the maximum eigenvalues of \H\, Eo, to the 
Coulomb corrected experimental energies [I3 of the 
lowest isovector-paired 0"*" states in even- A nuclei |33j| . 
Although the fits yield quite good agreement with the 
relevant experimental values, a reproduction of the fine 
properties of nuclear structure (typically of an order of 
magnitude or two less than the energies that were fit) 
is not guaranteed due to the strong mean-field contribu- 
tion. For the present investigation these parameters in 
the energy operator are not varied; their values are 
fixed as: G = 0.53rj, F = OAbil, C = 0.47, D = -0.97, 
E = -1.12(217), e = 9.36 in MeV for the I/7/2 level (with 
Ca core) and G = 0.3517, F = 0.3017, C = 0.19, 
-0.80, E = -0.49(217), e = 9.57 in MeV for the 

e) El. In the 



40 



D 

l/('5-)2p(| 3-)l5(9) shell (with a ^^Ni 
second case, the parameters of the effective interaction in 
the Sp{4:) model with degenerate multi-j levels are likely 
to be infiuenced by the nondegeneracy of the orbits. Nev- 
ertheless, as the dynamical symmetry properties of the 
two-body interaction in nuclei from this region are not 
lost, the model remains a good multi-j' approximation 
|33l | and the extent to which it provides for a realistic 
description can be further tested with the use of various 
discrete derivatives of the energy function. 



III. DISCRETE DERIVATIVES AND FINE 
STRUCTURE EFFECTS 

The symplectic 5*^(4) model [namely, the Eq maximal 
eigenvalues of \H\ reproduces the Coulomb corrected 
energies of the isovector-paired 0+ states quite well . 
A more detailed investigation and a significant test for 
the theory is achieved through the discrete approxima- 
tion of the d"^ Eq / dx"^ derivatives of the Eq energy func- 
tion. 



(m-l) 



m > 2, 



-J , m-even /^.n 



expressed recursively with two terminating conditions de- 
pending on the order m (even or odd) of the derivative, 
where the variable is a; = {n,i, N^i, N_i} according to 
the Sp{i) classification and (5 > 1 is a discrete integer 



step. The present investigation is focused predominantly 
on the (5 = 1 or 2 cases [the way the different variables, 
n,i,N+i or A^-i, vary in the 5*^(4) systematics can be 
recognized in Table ITT] . 

The first (m = 1) discrete derivative defined in jSJ, 
{Eo{x + 5) — Eq{x))/S, is related to the (5-particle sepa- 
ration energy, when x counts the (total, proton or neu- 
tron) number of particles. The general Stg'f^^ [x) quan- 
tity represents a finite difference between the Eq energies 
of neighboring nuclei, for example 

Stgf^ [X) = Eoi.+S)-2Eoi.) + Eoi.^S) ^ 

when m = 2, and 

(3)/ 



Stgr'ix) 



Ea{x+25)-3Eo{x+S)+3Ea{x)-Eo{x-5) /-x 
P ' 

when m = 3, and filters out contributions to Eq propor- 
tional to x"^"^ . 

The filters 10) are (m -I- l)-point expressions that ac- 
count for deviations from the common behavior of neigh- 
boring nuclei. When m > 3 the Stgg^\x) discrete 
derivative is independent of strong mean-field effects, 
strictly speaking it cancels out all regularly-varying lin- 
ear and quadratic in x contributions to the energy, that 
arc typically large, and only can provide for a descrip- 
tion of higher-order terms in the variable well as 
for discontinuities in the energy function. In this way, 
the finite energy difference isolates specific parts of the 
interaction that are comparatively smaller and may vary 
substantially from one nucleus to its neighbors. While 
these interactions do not contribute much to the overall 
trend of the Eq energies, they play a very significant role 
in determining nuclear structure properties. 

The mixed derivatives also provide useful information 
about the nuclear fine structure effects and are defined 
as 



Eo{x+Si,y+52)-Eoix+5i,y)-Eo{x,y+S2)+Eo{x,y) 
S1S2 



(8) 



where the variables represent quantities among the set 
{x,y) = {n, i, A^+i, iV_i} and 61^2 > 1 is a discrete incre- 
ment in accordance with the S'p(4) classification scheme 
(Table im. 

Different types of discrete derivatives are considered 
and various staggering patterns are investigated in the 
following sections. The corresponding components of the 
interaction isolated through the energy difference filters 
can be explained in analogous ways as in |23, ll^, in 
addition to the advantage that because they are free of 
Coulomb effect they reflect phenomena related only to 
nuclear forces. 



A. Discrete derivatives with respect to A'^+i and 
N-i: the N = Z region 

For even-even nuclei, the discrete approximation of the 
OEq /dN±i first derivative of the binding energies (in- 
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eluding the Coulomb repulsion energy) is related to the 
well-known two-proton (two-neutron) separation energy, 
which is usually defined as S'2p(2ri)(-^±i) = Eq{N±i) — 
E^{N±i — 2) [see Figure ^a) for a relation to proton 
number and Figure^b) for the difference of the Coulomb 
corrected energies Eq versus neutron number] . The <S'p(4) 
theor y r eproduces very well the available experimental 
data (33 (shown as 'x' or symbols for even-even nu- 
clei and as ' X- ' for odd-odd nuclei in Figure^, espe- 
cially the irregularity at A^+i = iV_i. The zero point of 
S2p along an isotone sequence determines the two-proton- 
drip line (dashed black line in Figure^, which according 
to the Sp{A) model for the l/(B-)2p(i a-^lgf^o^^ major shell 
lies near the following even-even nuclei: 



60 



Ge28, ^^Sego, ^«Kr 32, Sr34, Zr36, 
™Zr38, «2MO40, ''RU42, ^°Pd44, ^^Cd46, 



(9) 



beyond which the higher-Z isotones are unstable with 
respect to diproton emissions. These nuclei are not yet 
explored as seen in Figure^and an experimental compar- 
ison for the two-proton-drip line is expected to be soon 
possible due to radioactive beam experiments near the 
limits of stability. Yet, the findings of our model arc in 
close agreement with the results of other theoretical pre- 
dictions [H, \^ lil, li^. Particularly, the estimate for 
the two-proton separation energies in 0, con- 
firms the division in nuclides such that the isotones with 
lower/higher Z values than the nuclei in Q have pos- 
itive/negative S2p energies (compare to Figure^. In 
addition, the two-proton separation energies for those of 
the nuclei in ^ considered also in the other studies arc 
close in their estimates: the quadratic mean of the dif- 
ference in S2p between our model and ^3 is 0.32 MeV 
(in a comparison of the first three nuclei in (jni))jjs 0.78 
MeV when all the nuclei in ^ are compared to |4J| and 
is 0.43 MeV in a comparison to of the first four nuclei 
in Q. For odd-odd nuclei the zero point of S2p can be 
also determined (^OCasg, ^^^sai, ^Sfirsa, ^^Rb35, ^^^"37, 
78y39, ^"^Nbii, ^'^Tcis, ^°Rhi5, ^'^Agij) although it does 
not define the drip line, as S2p is a relation of the lowest 
0+ state energies Eq rather than of the binding energies 
for most odd-odd nuclei. 

As a whole, the higher-order derivatives with respect to 
proton (neutron) number have a smooth behavior. This 
is because these derivatives reflect changes only within 
a sequence of either even-even or odd-odd nuclei. The 
discretization of the d'^Eo/dN^.^ second-order derivative, 
4Wpp(„„)(Af±i) = Eo{N±i+2)-2Eo{N±i) + Eo{N±i-2) 

[= AStg^\N±i) 10 ], accounts for the interaction be- 
tween the last two pp {nn) pairs in the {N±i + 2) nucleus 
(Figure |2Ia)). The average interaction 6Ipp(^nn) may be 
used as an alternative way to 1281 to estimate the non- 
pairing like-particle interaction |50| [of the last two pro- 
tons (neutrons)]. It shows no outlined staggering pattern 
but a repulsive peak around the N = Z nuclei in very 
good agreement with the experiment |3^ and with the 
results and discussions of |28| . Another smaller peak is 
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FIG. 1: (Color online) The S2p two-proton separation energy 
in MeV for the even-even (ee) and odd-odd (00) nuclei in the 
l.f(5)2p(i 3)lg(9) rnajor shell: (a) versus number of protons 
for different isotones (A*" = 28 — 50) (the Coulomb repulsion 
energy is taken into account); (b) versus number of neutrons 
for Ge, Se, Kr, Sr isotopes. 



observed around midshell (Figure |2I a)), which is due to 
the particle-hole discontinuity introduced in the pairing 
theory. The analysis yields that as a whole the 5*^(4) 
model reproduces the fine structure effects in interactions 
isolated via the Stg'^ '{N±i) filters. 
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FIG. 2: (Color online) Second discrete derivatives of the Eq 
energy (l/^5j2pji ijlg^g-, shell): (a) with respect to A''+i, 
SIpp{N+-i_), as an estimation of the nonpairing like-particle 
nuclear interaction in MeV for the A'^ = 34, 36, 38 multiplets; 
(b) with respect to 7V+i and V_i, 5Vpn{N+i, N-i), for Zn, 
Ge, Sr isotopes. 

Another aspect of the nuclear interaction is revealed by 
the second-order discrete mixed derivative of the energy 
|4|,5T/p„(iV+i,iV_i) = (So(V+i+2,V_i+2)-£;o(iV+i + 
2,N-i)-EoiN+i,N-i + 2) + Eo{N+i,N-i))/4: ©. For 
even-even nuclei it was found to represent the residual 
interaction between the last proton and the last neutron 
[2^ 23 and it was empirically approximated by 40 /A 
|29| . The theoretical discrete derivative (Fi gur e |2Ib)) 
agrees remarkably well with the experiment [Sg, espe- 
cially in reproducing the typical behavior at A^+i = iV_i, 
and is consistent with the empirical trend (on average , 
~ 0.71 for 1/7/2 and ~ 0.52 for the major shell above the 
^^Ni core). It is well-known that the attractive peak in 
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the self-conjugate nuclei cannot be described by a model 
with an isovector interaction only |47j and in this respect 
our model achieves this result due to the additional terms 
included in the Hamiltonian, mainly the symmetry term 
(Figure (SJ- The 6Vpn energy difference provides for a 
powerful test for the symplectic model: the theory not 
only gives a thorough description of the isovector pn and 
like-particle pairing but additionally accounts for J > 
components of the pn interaction in a consistent way with 
the experiment. As a result the model can be used to 
provide for a reasonable prediction of SVpn of proton-rich 
exotic nuclei as well as odd-odd nuclei. 
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FIG. 3: (Color online) 5V^„ in MeV of the total binding en- 
ergy (■) and of the T = \ pairing energy (•) in comparison 
to experiment (x) for Ti-isotopes in the I/7/2 shell. The 
isovector pairing interaction is not enough to reproduce the 
experimental peak sk, N = Z . 



B. Discrete derivatives with respect to n and i: 
staggering behavior 

The <S'p(4) classification scheme can also be used to 
investigate energy differences with respect to the total 
number of particles n and their isospin projection i. In- 
deed, in contrast with the typical smooth behavior ob- 
served for discrete derivatives with respect to A^+i and 
iV_i that was highlighted in the preceding section, the 
derivatives with respect to n and i are the ones that 
reveal distinct staggering effects. They give a relation 
between even-even (ee) and odd-odd (00) nuclei and the 
patterns can be referred as an "ee — 00" staggering. 



1. Second and higher- order derivatives in one variable 

The discrete derivatives, Stg''i^\i), m = 1,2, show a 
prominent Ai = 1 staggering of the experimental energies 
|38l | of the lowest 0"*" isovector-paired states for different 
isobaric multiplets [see Figure 0] for the I/7/2 shell and 
Figure |3a) for nuclei above the ^^Ni core]. The theory 
reproduces this staggering very well. For each of the 
i multiplets (i fixed), a An = 2 stag geri ng effect is also 
observed for the experimental values [38| via the energy 
filters Stg2^^\n), m = 1,2, and successfully predicted 
by the symplectic model (Figure |H| (I/7/2) and Figure 
Uh) (l/(|)2p(i^3)lg(|))). 
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FIG. 4: (Color online) The Stg[^''^\i) discrete derivatives for 
different isobaric multiplets for even-A nuclei with valence 
nucleons in the I/7/2 shell with a core *°Ca. 





:(b) 

60 65 70 75 



85 90 95, 
A 



FIG. 5: (Color online) Discrete derivatives Stg^ {i) 
(l/(5j2pji i^lg^Qj, a '^^ Ni core): (a) 5 = 1, m = 2,3,4 for 
A = 76 isobars; (b) 5 = 2, m = 2,4 for (i = —1) multiplet 
[N = Z + 2\. 



The staggering amplitudes of both Stg'(^\i) and 
Stg^\n), while almost independent of the total num- 
ber of particles n, increase with increasing difference in 
proton and neutron numbers, i, and hence the ee — 00 
staggering effect is greater for the proton- (neutron-) 
rich nuclei than around N Z. Also, the amplitude 
of Stg^\i) increases in higher-order derivatives. This 
analysis shows a more complicated dependence of the en- 
ergy function on the isospin projection i than on the mass 
number A. 



The first. 



discrete derivative, S^, 



2Stg^i\n) 



m = 1, 

Eo{ji + 2) — Eo{n), where i is fixed, cor- 
responds to the energy gained when a, T = 1 pn pair is 
added [Figure E^a) (I/7/2) and Figure □(a ^^Ni core)]. 
Spn is the true pn separation energy only when Eq is 
the binding energy of the odd-odd nucleus involved in its 
calculation. The experimental data, where available |38l | , 
is also shown in Figure [7\ and the Sp{4) model follows 
the distinctive zigzag pattern very well. A An ~ 4 bi- 
furcation separates the nuclei into two groups: one of 
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FIG. 6: (Color online) Discrete derivatives Stgi^'^\n) for dif- 
ferent i multiplets for even-yl nuclei (I/7/2, a "Ca core). 



ferent isospin values (symmetry term). As noted in 
[27I l28l | , the significance of the various energy filters can 
be understood using phenomenological arguments that 
can be given a simple and useful graphical representa- 
tion. Specifically, each nucleus can be represented by an 
inactive core, schematically illustrated by a box, □, in 
which the interaction between the constituent particles 
docs not change. Active particles beyond this core can 
be represented by solid or empty dots, for protons or 
neutrons, above the box. 
The second-order filter 

Stg[^\i) = Eo{i + l)-2Eoii) + Eoii-l), 71 = const 
= Eo{N+i + l,N^i-l)-2Eo{N+i,N_i) 
+ So(iV+i-l,7V_i + l) (10) 



ses) 

|51j . Furthermore, the 
energy difference shows 



even-even nuclei [{n/2 + i)-cvcn] and another of odd- 
odd nuclei [{n/2 + i)-odd]. The Spn energy difference 
has a smooth behavior within each group. The mag- 
nitude of Spn is proportional to the total number of 
particles and increases (decreases) with i for odd-odd 
(even-even) nuclei (Figure O 
Stg^^\n) = ■S'*92"("+2)+g«g^"(»: 
no An = 4 staggering (average values of two consecu- 
tive data points in Figure [7^ . This indicates that the 
addition of an a-like cluster has almost the same ef- 
fect for both even-even and odd-odd nuclei. This state- 
ment does not contradict the stronger binding of even- 
pairs nuclei as compared to odd-pairs ones, which is 
detected via Spn and the binding energy (BE) filter, 

BEjZ + 2, TV + 2) - BE(Z+2,N)+BEiZ.N+2) ^ 



when centered at an odd-odd [{n/2 + i)-odd] self- 
conjugate {i = 0) nucleus, represents the pairing gap 
relation 2A 

(o\ , , m-~-odd) 00 •o 

stgY\i^o) = ' n + a-2a 

« 2A = 2App + 2A„„ - 4Ap„. (11) 

The result (|11() follows from the well-known definition of 
the empirical like-particle pairing gap 



= ^{BE{N+i ± 1, 7V_i T 1) - BEiN+i - 1, iV_i - 1) 
-2'[BE{N±i,N^i - l)-BE{N+i - l,7V_i - 1)])(12) 

= i("-n-2p -□]), 




FIG. 7: (Color online) The Stg'^\n) discrete approximation 
of the first derivative, dEo/dn {^^Ni core) with respect to: 
(a) A for several i multiplets; (b) i for different isobars. 



2. Pairing gaps 

The Stg["^\i) and Stg^\n) energy differences, m = 
1, 2, described above, isolate effects related to the var- 
ious types of pairing in addition to nonmonopole inter- 
actions resulting in changes in energy due to the dif- 



which isolates the isovector pairing interaction of the 
(iV±i)*'' and {N±i + 1)*'' protons (neutrons) for an even- 
even (iV+i — 1, A^_i — l)-core (marked by a square) [28l |. 
We also define the pn isovector pairing gap 



A„ 



= -{E„{N+,,N_, 



-[BE{K 
1 



BE{N+uN-i~i) 
+1 - 1, iV_i) - BEiN+i - 1, iV_i - 1)]) 

(13) 



= -(□-□-[□-□ 



as the pairing interaction of the (A^+i)*'' proton and the 
(A^-i)*'' neutron. In order to account correctly for the 
T = 1 mode of the pn pairing one should consider in H13|) 
the Eq energy of the odd-odd (Af+i,iV_i) nucleus (that 
is, the energy of the isobaric analog state rather than its 
ground state energy, BE). For the remaining even-even 
nuclei in (|l()|l . replacing the symbol Eq with BE is justi- 
fied. In the computation of A, all odd-A binding energies 
in (|12() and H13() cancel so their theoretical calculation is 
not required. 

The A relation of the gaps is a measure of the differ- 
ence in the isovector pairing energy between even-even 
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and odd-odd nuclei. For odd-odd N = Z nuclei infer- 
mation about A is extracted via the Stgl^'ii) energy fiL 
tcr 
A ? 



IT0|) . Both experimental and model estimations yield 
for all the odd-odd i = nuclei in the I/7/2 shell 



(for example, see solid (purple) line with open squares in 
Figure IHl for A — 46, i ^ 0). The result reflects the fact 
that in this case all three isovector pairing gaps, App, 



Ann and Ap„ , are equal [3l|, |3^ 
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FIG. 8: (Color online) Theoretical staggering amplitudes for 
the total energy in comparison to experiment |3g], for the 
isovector pairing energy, the pn and the like-particle pairing 
energies, and for the symmetry energy for A = 48, A = 46 
and A = 44 nuclei in the I/7/2 shell (a ''"Ca core). 

A different scenario regarding two aspects is encoun- 
tered when one considers the Stg\ (z) discrete derivative 
centered at an even-(n/2 + i) N = Z nucleus [relative to 
a (iV+i - 2,iV_i - 2)-core]: 



□ + □ -2 □ 



(- +i) even, (14) 



where an additional nonpairing two-body interaction 
jj^o,T=ii filtered out in this case. Here, for ex- 

ample, j^^^'^^^ is related to the nonpairing interaction 
of the three protons and of the three neutrons in the 

odd-odd nuclei H14|l . Another new feature of (|14|l is that 
(2) 

Stgi {i = 0) does not simply account for the energy 
gained when two pn pairs are created (in the first two 
odd-odd nuclei) and the energy lost to destroy a pp pair 
and a nn pair in the even-even N = Z nucleus. The 
straightforward reason is that pp, nn and pnT — \ pairs 
coexist. A good approximation that serves well in esti- 
mating the pairing gaps is to assume that a 2p — 2n for- 
mation above the inactive core (□) consists of no = 2/3 
pn pairs, ni = 2/3 pp pairs and n_i = 2/3 nn pairs 
[rather than a proton pair (ni = 1) and a neutron pair 
(n_i = 1)]. This is in analogy to an even-even n = 4 nu- 
cleus where the pp, nn and pn "numbers of pairs" are 
the same and equal to one third the total number of 
pairs, n/2 |3j|33. Additionally, the relations like Hll|) 
- H14|l are based on the assumptions that the interaction 



of a particle with the core is independent of the type of 
added/removed particles and is the same for all protons 
(neutrons) above the core. Finally, all the approxima- 
tions are of an order 0(l/ri). 

The additional nonmonopole two-body residual inter- 
action should be also taken into account for the 
rest i 7^ of the (ee, and 00) nuclei 



Stgf\i^Q) 



iA-f/,^^™ 



, ee 
, 00. 



(15) 



The main contribution to the I2 interaction is due 

to the symmetry energy as is apparent from the Sp{4) 
model. 

The very close theoretical reproduction of the ex- 
perimental staggering allows us to use the symplectic 
model as a microscopic explanation of the observed ef- 
fects through the investigation of the different terms in 
the Hamiltonian Q (Figure |HJ). According to the 5*^(4) 
model, the ee — 00 staggering patterns appear due to 
the discontinuous change of the seniority numbers driven 
by the T = 1 pairing interaction |33j . Even values of 
the seniority quantum number (z/i) in even-even nuclei 
and odd values for odd-odd nuclei lead to a change in 
pn and like-particle pairing energies in opposite direc- 
tions. After the contribution from the isovector pairing 
energy is taken away, the theoretical staggering ampli- 
tude, {—)^^^^^Stg^\i), has still a (typically large) com- 
ponent from the remaining {J ^ 0,T 1) interactions 
in the Hamiltonian Q), mainly the symmetry (T^) term 
[Figure |S1 long-dashed (purple) line with squares]. This 
is the same nonmonopole nuclear interaction, J^'^'^'"^^^, 
that was suggested in (|14|l and (|15|l using phenomenolog- 
ical arguments. Indeed, the symmetry energy contribu- 
tion is significant and nonzero in all nuclei but the odd- 
odd N = Z ^ (Figure IHl), which is consistent with the 
discussion above [ (|ll|l . (|14|) . H15|) ]. Also, an estimation 
of the pairing gaps is possible based on the examination 
of the model Hamiltonian but the theoretical staggering 
amplitudes of the T = 1 pairing energies (shown in Fig- 
ure IS)) need to be rcscaled in accordance with (|ll|l . H14|) 
and ifTCjl . 

In a way analogous to that used in (|15|l . the second- 
order discrete derivative with respect to n (can be com- 
pared to the filter used in [3ll |) 

(n) = Eoin+2)-2Eo^(n)+Eo(n-2) ^ i = COUSt (16) 

is related to the pairing gap relation 



Stg^^'in) 




,00, 



(17) 



where in the odd-odd case, for example, j^^'^-^^^ jg the 
non-pair interaction of the last two protons with the last 
two neutrons in the (n -I- 2) nucleus. The effects due to 
A cannot be isolated via (|17|l because of the additional 
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nonzero contribution due to the symmetry energy. How- 
ever, the staggering amphtude of the discrete derivative 
(|T?)|l . —3{—)^~^^Stg^\n), of the theoretical total, (nn) 
and pn pairing energies can provide for estimation of the 
pairing gaps A, App(„„) and — 2Aj,„, respectively [Figure 
IH^a)]. The like-particle pairing gap can be compared to 



the empirical value of A^ 



24/ A 



1/2 



ipai 

Q 



[sohd 



(purple) line]. The gap is smaller in odd-odd nuclei as 
compared to their even-even neighbors. This is a conse- 
quence of a decrease in the like-particle pairing energy in 
the odd-odd nuclei due to the blocking effect while there 
is an increase in energy due to the pn pairing. The pn 
isovector pairing gap increases toward i = and even- 
tually gets almost equal to App(„„) for odd-odd nuclei 
around the N = Z region, which is in agreement with 
the discussion of Isil l33l . 
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FIG. 9: (Color online) Estimation of the pairing gaps: (a) to- 
tal isovector pairing gap A, 2Ap„ and App -f A„„ , as well as 
the empirical like-particle pairing gap App -I- A„„ = 
shown for comparison, for ^4 = 48 and A = 46 nuclei 
versus i {If 7/2 shell); (b) like-particle pairing gap (accord- 
ing to 11811 ') versus A for i — ±6, ±7, ±8 multiplets in the 
l/(|)2p(i |)lsf(|) shell 



very large interaction matrices, as a very good approxi- 
mation one may use S'sym.T = ^^(^ + 1) with isospin 
values T = \i\ for even-even nuclei and T — |i|-t-l for odd- 
odd nuclei. Once the fourth-order discrete derivative 
of the approximated symmetry energy is removed from 
Stg'i'\i) fn^ . the like-particle pairing gaps App-|-A„„ are 
found to be in a very good agreement with the experimen- 
tal approximation of 24:/^/A for the {i = ±6, ±7, ±8) 
multiplets in the l/(5-)2p(i 3^1g(^a-^ major shell (Figure 
ini^b)). For lower |i| values the difference increases due to 
an increase in the pn pairing gap as mentioned above. As 
a whole, the agreement would not be possible if the sig- 
nificant energy contribution due to the symmetry energy 
was not taken into account. 



3. Second-order mixed derivatives 

Next we consider the second-order discrete mixed 
derivative of the relevant energies with respect to the 
total number n and the third projection i 



Stg'ilin.i) ^ 

Eo{n+2,i+l)-Ea{n+2,i)~Eo(nA+l)+Eo(n,i} 
2 



-|A + /,^^°'^^i 



, ee 
, 00, 



(20) 
(21) 



where in addition to the pairing gaps relation. A, there 
is the contribution due to the nonpairing interaction, 
I^'^^''^^^. For example, for the odd-odd {even-even} case 
it is the positive {negative} nonpairing average interac- 
tion between the last three protons {neutrons} in the {n+ 
2 {n},i + l) nucleus with a (n — 2 {n — 4},i) core. Within 
the 5*^(4) framework the additional nonpairing contribu- 
tion corresponds to the staggering of the symmetry en- 



ergy approximation. Eg 



of (- 



)"/2+.+i^(2|,| 



3). 



Furthermore, an average of the additional non-pair in- 
teraction is achieved by the fourth-order derivatives both 

-,(4)/ 



[Stg^^>{n)] and i 



3 

16 

3(- 



Stg[''\^)] 
{-T''+\Stg'f'ir) 



r/'+^{Stgi'\n) 



(18) 



//^°'^^i). (19) 



Assuming that the pn pairing gap is negligible for high- 
i nuclei in large shells like the l/(5j2p(i a^lf^j-g-) major 
shell, the gap relation ((TH|l or ifT^ll provides for a rough 
estimation of the like-particle pairing gaps. With the use 
of the model Hamiltonian we can estimate the addi- 
tional J^^^'^^^ interaction with the major input being 
the symmetry energy. Although the existence of a very 
small mixing of isospin values complicates the compu- 
tation of the symmetry energy for nuclear systems with 
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FIG. 10: (Color online) Second-order energy filter Stgl'l (n, i) 
for nuclei above the ^'^Ni core with respect to A (a) and i (b). 

The filter H2U|) isolates fine structure effects between 
two i multiplets [Figure ITUT a')] and two consecutive iso- 
baric sequences [Figure llOf b)]. Clearly, it reveals a 
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{An, Ai} = {2, 1} symmetric oscillating pattern as it is 
observed in the experiment |38l |. Its positive (negative) 
value is centered at even-even (odd-odd) nuclei and its 
amplitude increases (decreases) with |z|. This mixed dis- 
crete derivative (|20|l serves as another test for the 5*^(4) 
model and allows for a detailed investigation of the non- 
pairing, like-particle interactions involved. 

To isolate the effect of nonpairing interactions (again, 
it is understood to order 1/i^), an energy difference with 
respect to both N±i and i can be considered. The second 
discrete derivative of the energy 



z) = 

Eo{N±iA+l) + Eo{N±i,i), 



\22) 



Eo{N±i + lA+l)-Eo(N±i + l,i}- 
2 



represents the negative {positive} nonpairing two-body 
interaction of the last two neutrons {protons} with a pro- 
ton and a neutron in the {N±i + l,i {+1}) nucleus. It 
shows prominent Ai = 1 staggering patterns for different 
i multiplcts fFigurc [TT|l . While in the framework of the 
Sp{A) model its amplitude does not depend on N±i and 
i except for irregularities around the mid-shell, the mag- 
nitude of the few experimental values [s^ (where data 
exist) tends to be slightly lower away from the closed 
shell. As a whole, the results show that the staggering 
behavior of this interaction is due to the fine structure 
features in the relationship between the like-particle and 
pn nonpairing interactions and differs between proton- 
rich and neutron-rich nuclei. 
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FIG. 11: (Color online) Discrete derivative, Stgfl{x,i), for 
various i multiplets for even-^ nuclei; (a) x = A'^+i, I/7/2 
level; (b) x = A''+i, l/(5-)2p(i _3)1(?(9) shell; (c) x = N-i, 
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Regarding H22() and the other discrete approximations 
of the derivatives in Section Fill Bl it is clear that the os- 
cillating patterns that exist and their regular appearance 
throughout the nuclear chart cannot be a simple artifact 



due to errors in the experimental or theoretical energies. 
Even more, the staggering amplitudes are usually (very) 
large compared to the energy uncertainties. 

For all the discrete derivatives that we have investi- 
gated above and that show "ee-00" staggering behavior, 
the discontinuity of the symmetry term (due to discrete 
changes in the isospin value) plays an important role. In 
contrast, when these discrete derivatives include states 
of odd-odd nuclei with a dominant T = pn coupling 
there is a constant or no contribution due to the symme- 
try energy, and hence yield patterns of different shapes 
and interpretations. Our investigation does not aim to 
account for such effects. It is focused on the ee — 00 stag- 
gering behavior of the £^0 energies of the lowest isovector- 
paired states as observed from the experimental data and 
reproduced remarkably well by the Sp{A) model. 



IV. CONCLUSIONS 

A dynamical 5*^(4) symmetry was used to provide for 
a natural classification scheme of nuclei and to describe 
isovector pairing correlations and high- J interactions. In 
a previous study [23| , it was found that the 5*^(4) model 
reproduced reasonably well the experimental energies of 
the lowest isovector-paired 0^ states and provided for an 
estimation of the interaction strength parameters. 

Here the sp(4) algebraic approach has been further 
tested through second- and higher-order discrete deriva- 
tives of the energies of the lowest isovector-paired 0"*" 
states in the 5'p(4) systematics, without any parameter 
variation. If reality were only a mean-field theory, none of 
the finite energy differences would reveal regular or irreg- 
ular staggering effects. The reason is that any effect due 
to a smoothly varying mean-field part of the nuclear in- 
teraction is either entirely cancelled out in a finite energy 
difference filter or contributes regularly to the isolated 
part of the interaction. Indeed, the results obtained show 
that this is not the case and staggering behavior is ob- 
served. The theoretical discrete derivatives investigated 
not only followed the experimental patterns but their 
magnitude was found to be in a remarkable agreement 
with the data. The proposed model successfully inter- 
preted the following: the two-proton (two-neutron) sep- 
aration energy S2p{2n) (hence determined the two-proton 
drip line) for even-even nuclei, the Spn energy difference 
when apnT = 1 pair is added, the observed irregularities 
around N ~ Z, the prominent ee — 00 staggering when 
even-even and odd-odd nuclides are considered simulta- 
neously, the like-particle and pn isovector pairing gaps, 
and the large contribution to the finite energy differences 
due to the symmetry term. The oscillating effects, where 
observed, were found to develop due to the discontinu- 
ity of the seniority numbers for the pn and like-particle 
isovector pairing, which is in addition to the larger stag- 
gering due to the discontinuous change in isospin values 
(symmetry term) between even-even and odd-odd nuclei. 

We found a finite energy difference that, for a spe- 
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cific case, can be interpreted as an isovector pairing gap, 
A = App + Ann — 2Ap„, which is related to the Uke- 
particle and pn isovector pairing gaps. They correspond 
to the T ~ 1 pairing mode because we do not consider 
the binding energies for all the nuclei but the respective 
isobaric analog 0+ states for the odd-odd nuclei with a 
J ^ 0+ ground state. This investigation is the first of 
its kind. Moreover, the relevant energies are corrected 
for the Coulomb interaction and therefore the isolated 
effects reflect solely the nature of the nuclear interaction. 

The outcome of this investigation shows that, in com- 
parison to the experiment, the simple 5*^(4) model repro- 
duces not only global trends of the relevant energies but 
as well the smaller fine structure effects driven by isovec- 
tor pairing correlations and higher- J pn and like-particle 
nuclear interactions. In particular, the sp{A) algebraic 
model was used to interpret specific phenomena revealed 
in finite energy differences and to investigate the con- 
tribution of the underlying interactions. In this way, it 
provides for an estimation of the isovector pairing gaps. 
For N = Z odd-odd nuclei all three pairing gaps were 
found equal while the pn pairing was found to weaken rel- 
ative to the likc-particlc pairing strengths with increasing 



proton (neutron) excess. The like-particle pairing gaps 
were found to be in a good agreement with the empiri- 
cal value of 12 /^/A. Additionally, the discrete derivatives 
give insight into particular small parts of the various non- 
(J = 0,T = 1) interactions, mainly into the detailed 
contribution of the interaction related to the T{T + 1) 
term (symmetry energy). Small deviations from the ex- 
perimental data are attributed to other two-body inter- 
actions or higher-order correlations that are not included 
in the theoretical model. 

We explored independent finite energy differences 
based on a simple sp(4) algebraic classification scheme. 
The results suggest that this theoretical framework can 
be used to reproduce various experimental results includ- 
ing observed staggering behavior in fine structure effects 
of nuclear collective motion. 
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